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APPARATUS FOR SIMULATING THE DEFORMATION OF 
MATERIALS, NOTABLY OF SOFT BODY TISSUES 

CROSS-REFERENCE TO RELATED PATENT APPLICATIONS 

[0001] This application claims priority to French Patent Application 0101717, filed 
02/08/2001 including the specification, drawings, claims and abstract, incorporated 
herein by reference in its entirety. 

BACKGROUND OF THE INVENTION 

[0002] The following description of the background of the invention is provided 
simply as an aid in understanding the invention and is not admitted to describe or 
constitute prior art to the invention. 

[0003] The invention relates to the simulation of deformation of materials, notably 
of soft body tissues. 

[0004] Usually, a display interface is provided for representing at least one object 
following a predetermined shape, displayed with a selected grid pattern. A user 
interface provided with mechanical handling means allows the user to simulate one or 
more forces exerted globally on the object. 

[0005] A memory is also provided for storing mechanical parameters of the material 
of which the object in question consists, data as to the position of the object, recorded 
at the vertices of a selected mesh, and force data which represent, in intensity and 
position, stress to be exerted on the object. 

[0006] A computer operatively connected to the memory evaluates new positions of 
the vertices, as a function of the stress exerted and the mechanical parameters of the 
material. In order to display the deformation and/or the return to a resting position of 
the moving object, the computer calculates the evolution of the positions of the 
vertices over time as a function of this stress. When it is used to simulate the 
deformation (or cutting) of body tissues, a qualified person (such as a surgeon) is able 
to follow the simulation of a deformation and/or a return to a resting position. 
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[0007] In the published application no. FR-97 14506 by this Applicant, the 
simulation of the deformations of the object was based on a linear model, in the 
embodiment described. Although the use of the apparatus described in FR-9714506 
has proved very promising in the field of surgical assistance for simulating the cutting 
of soft body tissues (liver, skin, etc) , the simulation of large movements, possibly 
non-linear, was not completely natural. Moreover, this apparatus could not predict 
the reaction to a force exerted on just any type of deformable material, particularly in 
the case of large movements. 

[0008] Thus, simulated deformations with large shifts in position are not normally 
reliable. The present invention sets out to improve the situation. 

SUMMARY OF THE INVENTION 

[0009] According to one embodiment of the invention, the computer comprises a 
module capable of repeatedly determining, for each mesh element, the deviation 
between the current length of a side of the mesh and its resting length, and of storing 
the respective force data for each vertex of the mesh, and a module capable of 
repeatedly determining, for each vertex, new positional data relating to this vertex as a 
function of the composition of the forces exerted thereon, data relating to at least one 
earlier position of the vertex and mechanical parameters of the material. 

[0010] According to one embodiment of the invention, the mesh pattern chosen 
comprises at least one triangular base. 

[0011] Thus, if the object is in the form of a hollow three-dimensional envelope, the 
mesh pattern chosen is preferably triangular, so that the computer is capable of 
determining the composition of the forces at each vertex of a triangle, as a function of 
the deviation between the current length of each side of the triangle and the length of 
this side at rest. 

[0012] If the object has a solid three-dimensional shape, the mesh pattern chosen is 
preferably tetrahedral, so that the computer can determine the composition of the 
forces at each vertex of the tetrahedron, as a function of the deviation between the 
current length of each edge of the tetrahedron and the length of this edge at rest. 
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[0013] Advantageously, the memory is arranged so as to store, in conjunction with 
each mesh element, mechanical parameters of the material of the mesh, at least partly 
defined locally, particularly at the level of the mesh or elements thereof. 

[0014] According to another embodiment of the present invention, the computer is 
designed to determine said deviation between current length and resting length, in 
order to estimate a derivative of the potential energy of deformation of each mesh, 
expressed as a function of a Green-St Venant tensor and mechanical coefficients 
inherent in the material in this mesh. 

[0015] The derivative of the potential energy thus expressed in relation to the 
position of the vertex makes it possible to obtain the force (or composition of forces) 
exerted thereon. 

[0016] In estimating the derivatives of the potential energy at the vertices, thus 
expressed as a function of a Green-St Venant tensor, the computer is adapted to 
determine more precisely the differences between the squares of the current lengths 
and resting lengths of each side, in order to determine the composition of the forces 
mentioned above. 

[0017] The mechanical parameters mentioned above advantageously comprise the 
Lame coefficients of the material in each mesh in question. 

[0018] The apparatus according to one embodiment of the invention, can thus be 
used to simulate large movements and/or non-linear elasticity, possibly of 
heterogeneous objects, formed of different materials. 

[0019] It is thus possible to simulate the deformation of a solid three-dimensional 
object, possibly enveloped in a hollow skin. In this case it will be advisable to apply a 
tetrahedral grid to the object and a triangular grid to the skin and to apply suitable 
mechanical parameters to each type of mesh. It is also possible to simulate the 
movements of soft objects (such as muscles) in the presence of rigid integrated 
objects (such as bones), as will be seen hereinafter. 

[0020] The present invention may also take the form of a process for simulating the 
deformation of objects, particularly of soft body tissues, comprising all or some of the 
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following steps of (a) specifying a predetermined form for the object at rest, (b) 
making a grid image of said form, based on a selected pattern, preferably with a 
triangular base (triangle or tetrahedron), (c) storing the positional data for the vertices 
of the mesh in the object at rest, (d) simulating one or more forces exerted globally on 
the object (or the movements of the object), (e) estimating, as a function of the 
mechanical properties of the material and the positions of the preceding vertices, a 
resultant force on each vertex of each mesh, and (f) as a function of the force exerted 
on each vertex, determining the movement and, hence, the deformation and/or the 
return to a position of equilibrium of the vertices of the meshes. 

[0021] Step (c) also involves attributing indices to the elements of the meshes, 
within the scope of establishing a data structure, and, if appropriate, simplifying the 
structure of the object, as will be seen hereinafter. Steps (e) and (f) are based on a 
model developed by the Applicant which will be described in detail hereinafter. 

[0022] It is to be understood that both the foregoing general description and the 
following detailed description are exemplary and explanatory only, and are not 
restrictive of the invention as claimed. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0023] Further features and advantages of the invention will become apparent from 
studying the detailed description that follows and the accompanying drawings 
described below. 

[0024] Figure 1 diagrammatically shows an apparatus according to one embodiment 
of the invention. 

[0025] Figures 2A and 2B diagrammatically show triangular meshes, in the initial 
and deformed state, respectively, with notations for the points of the triangle 
subsequently used in estimating the composition of the forces at each vertex of the 
triangle. 
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[0026] Figure 3A is a diagram of a tetrahedron which helps in the understanding of 
the notations used for a point Pj of the tetrahedron and the normal vector nj of the 
opposite face. 

[0027] Figures 3B to 3D are diagrams of tetrahedrons which, together with the 
abovementioned Figure 3 A, help to show how the sums on the edges and pairs of 
edges of a tetrahedron are calculated, in order to evaluate a force to which a vertex Pp 
of the tetrahedron is subjected. 

[0028] Figures 4A and 4B are diagrams to assist in the understanding of the 
attribution of indices in different triangular meshes. 

[0029] Figures 5 A and 5B are diagrams of triangular motifs comprising, 
respectively, a current vertex (represented by a square) and a virtual vertex 
(represented by a disc), before simplification, and the same gridded shape after 
simplification. 

[0030] Figure 6 diagrammatically shows an example of a nonconforming triangular 
mesh. 

[0031] Figure 7 shows a diagrammatic flow chart of the progress of the calculations 
of forces and positions of the vertices representing the object. 

DETAILED DESCRIPTION 

[0032] The present description will refer to formulae and equations which appear in 
the annex. 

[0033] The drawings, and the annex, essentially contain elements of a definite 
nature. They can therefore not only be used to gain a fuller understanding of the 
description, but may also contribute to the definition of the invention, where 
appropriate. 

[0034] Reference is made first of all to Figure 1 in which the apparatus is shown in 
the embodiment described in the form of a computer (PC). The computer comprises a 
permanent memory MEM (for example of the ROM type) in which are stored 
programs for calculating the positions of the mesh vertices and the resultant forces at 
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these vertices, as well as a processor jliP operatively connected to a working memory 
MT (for example of the RAM type) for carrying out these calculations. 

[0035] The abovementioned modules of the computer described hereinbefore may, 
if necessary, be in the form of programs stored in the memory MEM, or may be 
loaded into the working memory MT, from reading a data carrier in an external 
peripheral (not shown) such as a disc drive or a CD-ROM reader, or received via a 
communications link, e.g. a MODEM. Of course, such programs appear to be 
particularly important means for carrying out the invention. Thus, the present 
invention also relates to programs of this kind, possibly on a data carrier of the type 
mentioned above (disc, CD-ROM, etc), detailed embodiments of which are described 
hereinafter. 

[0036] The computer further comprises a monitor VISU connected to the processor 

via a graphic interface IG, the unit consisting of VISU and IG (referred to as the 
"display interface" hereinbefore) providing a display of an object the deformation of 
which is being simulated. The computer further comprises a user interface IU 
connected to mechanical handling means such as a keyboard CLA or a mouse MO, 
enabling a user to simulate a force exerted globally on the object, or a change of 
position of this object, with initial deformation. 

[0037] The shape of the object to be displayed in grid form may be a surface shown 
in three dimensions. In this case, the mesh is triangular (figure 2 A). The shape of the 
object to be displayed in grid form may be a three-dimensional solid and in this case 
the mesh is tetrahedral (Figure 2B). If appropriate, a grid module is stored in the 
memory MEM and can be implemented by the processor [iP, co-operating with the 
graphic interface IG, to display a shape of the grid- covered object on the monitor 
VISU. 

[0038] The following is a description of a model, for example purposes only, 
developed by the Applicant and applied to an apparatus for simulating deformation. 
This model uses a resting form without mathematical induction in relation to a resting 
position (mathematical induction by shifts which is normally the basis of conventional 
models). 
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[0039] Let us consider first of all a surface shown in three dimensions, gridded with 
a triangular mesh pattern. A triangle TO = (P0, PI, P2) of the resting form, with edges 
(loJhh)- This triangle (on the left in figure 2A) undergoes a transformation <D to 
give the triangle T = (QO, Ql, Q2) (on the right in figure 2B). 

[0040] The potential elastic energy of this deformation is expressed here from the 
Cauchy-Green C deformation tensor (formula Fl of the Annex). Looking at the 
finished element obtained by associating a set of linear interpolation functions with 
this triangle, the deformation gradient is expressed at every point X on the triangle by 
the formula F2, wherein A is the area of the triangle TO, nj is the unit vector 
perpendicular to the segment (Pj+1, Pj+2), and the aj's are shape vectors of the 
triangle. 

[0041] The deformation tensor C is then given by the expression of formula F3. 
Advantageously, this deformation tensor is invariant by translation and therefore does 
not depend on the choice of the origin of the reference point. Preferably, the origin of 
the reference point is placed in the centre of the circle circumscribed around the 
deformed triangle T, of radius r (figure 2A) , which allows the tensor C to be 
expressed in the form given in formula F4, where Lp is the length of the deformed 
edge joining the vertices Qj and Qk. 

[0042] The expression of the potential deformation energy (given by formula F5 in 
the annex) uses those of the trace of the tensor C (given by formula F6) and of the 
trace of the tensor C2 (which will be described hereinafter). 

[0043] Starting from the classic equations in the triangle TO (cosine rule), which are 
expressed by the formula F7, the scalar products between the directing vectors are 
obtained, given by equation Rl in the annex. The respective traces of the tensors C 
and E are thus written according to formulae 8 and 9 in the annex. This term depends 
only on the variations in length of the edges of the triangle and its resting 
shape(represented by the magnitudes Ki). 

[0044] From the expression of tr C 2 (formula FIO), trE 2 is deduced and, finally, the 
potential elastic energy of deformation of the triangle is given by equation R2 in the 
annex. 
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[0045] The elastic force applied by the triangle at one of its vertices is obtained by 
deriving the elastic energy relative to the position of this vertex. In the embodiment 
shown in Figure 2A, the vertex P0 has QO for its deformed position. Deriving the 
equation R2 in relation to QO, the first step is preferably to derive the terms L p in 
order to obtain the expression of the force according to formula Fll, which is the 
equation of the elastic force exerted by the triangle at each of its vertices. This force 
is expressed solely as a function of the variations in length of the sides of the triangle 
(Li 2 - li 2 ) , the mechanical parameters of the material (X and \x) and the shape of the 
triangle at rest (given by the terms Ki). 

[0046] Advantageously, this force implicitly contains a property of rotationaly 
invariance as all its terms remain zero as long as the lengths of the sides of the 
triangle do not change, i.e. as long as the triangle is not deformed. Thus, the 
calculations of movements can be simplified by virtue of the fact that it is no longer 
necessary to compensate for rigid transformations, as invariance by rigid 
transformation is intrinsic to the formulation of the deformable model. 

[0047] Extending the foregoing to a solid three-dimensional shape (having a 
volume), gridded with a mesh of tetrahedrons, by formulating the elasticity according 
to Saint Venant-Kirchhoff it is possible to express the gradient of deformation 
(formula F12), as well as the Cauchy-Green deformation tensor (formula F 13) , where 
Aj is the area of the face opposite the vertex Pj, nj is the perpendicular to this face and 
V is the volume of the tetrahedron (Figure 3A). 

[0048] By placing the centre of the reference point at the centre of the sphere 
circumscribed around the tetrahedron T, it is possible to express the scalar product of 
the position vectors of two vertices of the tetrahedron as a function of the length of the 
edge that joins them and the radius r of the sphere, thus making it possible to write the 
deformation tensor according to formula F14. 

[0049] The elastic deformation energy is expressed as a function of the Green-Saint 
Venant deformation tensor (bearing the reference E), but in a 3D grid (formula F15). 
The expression of the trace of the tensor C is also expressed with the scalar products 
(aj.ctk), using the cosine rule for dihedral angles and as a function of the lengths of 
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the edges of the tetrahedron (formulae 16 and 17). Consequently, the trace of C 2 is 
written according to formula PI 8. Finally, the Applicant has shown that the elastic 
deformation energy is expressed by the formula F19 of the annex (complemented by 
the formula F20). 

[0050] The expression of the elastic force exerted by the tetrahedron on one of its 
vertices is obtained by deriving the elastic energy relative to the position of this 
vertex. A formulation which is easy to use is expressed by formula F21 of the annexe 
and is interpreted using Figures 3A to 3D. 

[0051] In formula F21 

in the second line, k and m are the complementary indices of j and p, 

in the last line, m is the complementary index of j, k and p, and : Kjkmp = 
KjpKkm + Kj m Kk p -KjkK mp 

[0052] The total on the edges (j,p) indicates that it is necessary to travel over the 
three edges adjacent to the vertex Pp (on which the force is exerted) whereas the total 
on the edges (j,k) (indices different from p) involves the three edges located on the 
face opposite the vertex Pp. 

Referring to Figures 3 A to 3D: 

- the total on the edges (j,p) involves three terms (Figure 3 A), 

- the total on the pairs of edges ((j,p), (k,p)) involves three terms (Figure 3B), 

- the total on the pairs of edges ((j 5 p), (j>k)) involves six terms (Figure 3C) and 

- the total on the pairs of edges ((j ? k)> (m,p)) involves three terms (Figure 3D). 

[0053] Once again, there is rotational invariance of the above expressions, for a 
solid three-dimensional form, which advantageously makes it possible to avoid the 
appearance of unrealistic distortions when part of the object is rotated too far. This 
invariance also extends to movements of translation. Thus, the object may undergo 
any overall movement of translation and rotation while retaining zero internal energy 
and stresses. The Applicant has verified that in this case the rigid component of the 
deformation is zero. 
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[0054] Moreover, the discretization of the object using finished elements comes to 
approximate the elastic energy over a set of sub-domains defined by simple geometric 
shapes. The above model thus suggests, for each sub-domain, a formulation of the 
energy which only depends on the change in shape of this domain. 

[0055] Advantageously, displaying a simulated deformation while considering an 
initial resting shape (rather than initial resting positions) appears to be more natural 
than the deformations simulated with known models. Moreover, considering an 
initial resting shape avoids the problems of accuracy that arise when the model moves 
substantially away from its resting position. 

[0056] The following is a description of the preliminaries to using the 
abovementioned computer modules to determine the forces at each vertex of the 
chosen mesh. 

[0057] First of all, a data structure should be defined. The Applicant has first dealt 
with the difficulty of attributing suitable indices to the elements of each mesh. 

[0058] For triangular patterns, the data structure corresponding to triangulation 
consists of a list of triangles, a list of sides and a list of vertices: 

- for each vertex, an adjacent side is determined, 

- for each side, a pointer is determined on the two adjacent triangles, and the two 
vertices which form the side are determined, 

- for each triangle, the pointers are determined on the three vertices that form the 
triangle, the three sides of the triangle are determined, the three triangles adjacent to 
each side are determined, and three values giving the index of adjacency of the 
current triangle for each adjacent triangle are determined. 

[0059] The orientation of the three vertices in the triangle is carried out as follows. 

[0060] For the grid as a whole to be correctly oriented, each edge or side must be 
traveled over in two different directions when each of its two adjacent triangles is 
traveled over. Referring to Figure 4A, the side Ai of index i is the side opposite the 
vertex Si of index i in this triangle. The adjacent triangle to index i is the neighboring 
triangle TV, adjacent to the side Ai. Finally, each triangle always has three adjacent 
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triangles. In order to represent surfaces with borders (a plane has one border, a 
cylinder has two borders) the notion of so-called "virtual" vertices and triangles is 
used (Figure 4B) . The virtual triangles are not displayed and the position of the 
virtual vertex is of no consequence to the geometry of the grid. 

[0061] An example of index attribution will now be given by way of example, in 
C++ language. 

[0062] If the adjacent triangle number i is called neighbor [i], then the adjacency 
index "whichNeighbor[i]" is defined by an instruction such as 

Neighbor[i] — >neighbor[whichNeighbor[i]]=this ; 

[0063] For reasons connected with a tetrahedral grid, the data structure of a 
tetrahedral volume consists of a list of tetrahedrons, a list of triangles, a list of edges 
and a list of vertices. 

[0064] For each vertex, its position and an adjacent tetrahedron are stored. 

[0065] For each edge, the pointers on the two adjacent vertices and a pointer on a 
triangle adjacent to this edge are stored. 

[0066] For each triangle, the pointers on the three vertices that form the triangle and 
the pointers on the two tetrahedrons adjacent to the triangle are stored. 

[0067] For each tetrahedron: 

- the pointers on the four vertices that form the tetrahedron, 

- the pointers on the six edges of the tetrahedron, 

- the pointers on the four triangles of the tetrahedron, 
and 

-the pointers on the four adjacent tetrahedrons opposite each of the four vertices are 
stored. 

[0068] Each tetrahedron is oriented. Starting from the orientation of the four 
vertices ABCD, the operator seeks to orient the four triangles : BCD, ADC, ABD, 
ACB. To ensure correct orientation of the grid, each triangle has to be traveled over 
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in two different directions by traveling over each of its two adjacent tetrahedrons. 
Finally, each tetrahedron is always adjacent to four tetrahedrons. The notion of so- 
called "virtual" vertices and "virtual tetrahedrons" is introduced to take account of the 
borders of the grid. 

[0069] In order to simplify the structure of the deformed object, the grid is cut up by 
removing material (more precisely a vertex Sr, although this vertex cannot be a virtual 
vertex). 

[0070] In the case of a triangular mesh, if this vertex is not adjacent to a virtual 
vertex, this vertex is regarded as being virtual like all these adjacent triangles. If the 
vertex is already adjacent to a virtual vertex, the side which connects this point to the 
virtual vertex is removed, and all the triangles adjacent to the virtual vertex are forced 
to be virtual (Figures 5A and 5B). 

[0071] With a tetrahedral grid pattern, the procedure involves removing a set of 
associated tetrahedrons. Once the tetrahedrons have been removed the grid must 
remain true, as described hereinafter. For this, it may be advantageous to remove 
more tetrahedrons than those present in the set of tetrahedrons mentioned at the 
beginning. Virtual vertices may also be added or removed when creating associated 
components or removing the borders. 

[0072] The Applicant has confronted the problem of ensuring conformity of the 
grid, a definition of which is provided hereinafter. 

[0073] For a triangular mesh to be true, each side must belong to more than two 
triangles and the topological neighborhood of each vertex must remain 
homeomorphous with a circle (or with a semicircle), i.e. if two triangles share the 
same vertex, there must be a set of triangles adjacent to each other which make it 
possible to go from one triangle to another. Referring to Figure 6, the triangular mesh 
shown is not true, as the two triangles have an intersection along a vertex, while there 
are no adjacent triangles making it possible to go from one of these two triangles to 
the other. 
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[0074] For a tetrahedral grid to be true, each triangle must belong to more than two 
tetrahedrons and the topological neighborhood of each vertex must be 
homeomorphous with a sphere (or a hemisphere if it is on the border). Therefore, if 
two tetrahedrons share the same vertex, it is essential that there is a set of tetrahedrons 
adjacent to one another which make it possible to pass from one tetrahedron to 
another. 

[0075] The equation that governs the movement of the linear elastic model is the 
Newtonian differential equation El in the annex. This equation is connected with the 
differential equations encountered in the mechanics of continuous media (equation 
E2). 

[0076] According to the theory of the finished elements, the mass matrix M and the 
attenuation matrix C are hollow and depend on the properties of each tetrahedron. 
Moreover, these matrices M and C may be considered here as being diagonal, i.e. the 
mass and the effects of attenuation are concentrated on the vertices of the mesh. This 
simplification, known as mass-lumping, uncouples the movements of each of the 
nodes and thus makes it possible to write the previous equation E2 as a set of 
independent differential equations of the type El, with an equation for each vertex. 

[0077] The choice of an explicit integration diagram leads to the expression of the 
elastic forces at the time t, to calculate the new position of the vertex at the time t+1, 
with the equation E3 in the annex. Advantageously, this equation E3 applies to all 
types of force, both linear and otherwise. 

[0078] The explicit diagrams are only conditionally stable and therefore converge 
less rapidly than the implicit diagrams. 

[0079] Another difficulty arises : for the diagram to be stable, the time step used for 
the digital integration should be sufficiently short. The critical time step depends on 
the greatest intrinsic value of the rigidity matrix and of the local values of the mass 
and of the attenuation. To optimize the time step, the mass matrix is fixed to be equal 
to identity, and the optimization of the local values of attenuation, so as to limit the 
risks of oscillations, have the effect of improving the packaging of the linear system, 
and hence of establishing its resolution. 
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[0080] Reference will now be made to Figure 7, wherein a flow chart 
diagrammatically shows the different treatment steps applied to a triangulated three- 
dimensional surface, by way of example. 

[0081] The treatment begins with a first step 10, in which the structure of the data 
described above is applied, with a suitable attribution of indices. Then in step 12 
consideration is given to each triangle T n , for which the Lame coefficients X n and yi n 
are sought, as well as the area of the triangle A (step 14). For each side of the triangle 
(j being an integer from 0 to 2), the difference between the squares of the lengths of 
the sides is calculated, between the current length lj and the previous length Lj of the 
side(step 16). From these differences in the lengths of the sides squared, the forces 
exerted at each vertex Sj of the triangle T n are deduced (force Fj n calculated in step 1 8 
of the flow chart). Moving on to the next triangles, the forces are integrated over all 
the triangles by looping on the index n (step 20). 

[0082] In steps 30 and 32, the new positions Q' n of the vertices Sj of the deformed 
triangles are determined, from their previous position Q n and the forces F n . By 
looping on to the index n in step 36, integration is carried out over all the deformed 
objects. The new positions of the vertices Sj are preferably estimated by using the 
equations El, E2 and E3 mentioned hereinbefore. 

[0083] In the flow chart in Figure 7, the two respective integrations (of the forces Fj n 
and of the positions of the vertices Sj in the deformed objects) are shown separately. 
Of course, it may be that some of the steps 12 to 20 may be carried out at the same 
time as or after some of the steps 30 to 36. This choice depends on the structure of 
the two modules described hereinbefore, contained in the abovementioned computer, 
and, if applicable, on the data structure. 

[0084] A test 42 is preferably provided to find out whether the object has fully 
returned to its resting position, i.e. : is there an edge in the object the length of which 
differs substantially from the previous one? 

[0085] If the computer does not find such an edge as it travels over the triangles, the 
object has returned to its initial resting shape and the treatment is stopped at 40. 
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[0086] Otherwise, an image of the object is displayed by the monitor VISU 
mentioned above and a time loop (step 44) is incremented, to pass to the next time 
and show a progressive return to a resting shape of the object, using equation E3, if 
appropriate. This time loop can be initiated after the first starting step 10. 

[0087] Of course, the invention is not limited to the embodiment described 
previously by way of example, but also extends to other alternative embodiments 
defined within the scope of the claims that follow. 

[0088] Thus, it will be understood that the invention is not limited to an application 
to simulating the deformation of soft body tissues. As described previously, the shift 
in position of a rigid object (for example a bone attached to a muscle) integral with 
the soft tissue can also be simulated. In some cases it may be appropriate to use a 
tensor (of inherently known expression) rather than a deformation tensor of the type 
described above. 

[0089] The formula for the deformation energy W, in both a triangular and a 
tetrahedral grid, is provided hereinbefore by way of example and depends on the 
material of which the object whose deformation is being simulated consists. 
However, it is generally expressed by a combination of the traces of the tensor E and 
the tensor E 2 . In these particular cases, only the trace of the tensor E (or the tensor E 2 ) 
may be used. 

[0090] The foregoing description of a preferred embodiment of the invention has 
been presented for purposes of illustration and description. It is not intended to be 
exhaustive or to limit the invention to the precise form disclosed, and modifications 
and variations are possible in light of the above teaching or may be acquired from 
practice of the invention. The embodiment was chosen and described in order to 
explain the principles of the invention and as a practical application to enable one 
skilled in the art to utilize the invention in various embodiments and with various 
modification are suited to the particular use contemplated. It is intended that the 
scope of the invention be defined by the claims appended hereto and their equivalents. 
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ANNEX. 



Fl : <7=V#*V$. 

i=0 j=o 2** 



F2 



F7 



F3 = (? = I]CQi-Qib)(o: i ®Q: jfc ) 



F4 : C = 



Er 2 («, *«,) + E ~ f ) ® «*> 



- -4e 2 ?(°*®°'*) 

F5 : W = ^(trJS^ + zrfr^ 2 

trJB 2 = \\tr C* \-trE-l 
41 1 2 



F6 *rC?=-^E X ?( a i-°*) 



3& 



{lj h co$($ P ) = \ {1} + H - 5) = K,, if S jt r , j = coroptp, (?) 
an&fc = eomp(p,r) 
V if ff = r , j = comp(p,q) 



Rl 



£42 iifj = fc 
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F8 



F9 



FIO 



Fll : 



9Qo 



F12 



F13 : 



F14 



+ M((^-/f)J2(l|Q^-^oQrQi)+(x|-^)/3(|f Q^ -KoQ^) 
-(^o - $ (if *i QiQ$ +lf A qbQ2)) J 



4* 



07 = X;(Qj. <8><* fc ) 



-{;-* 



if J = fc 



F15 



trB 
trE 2 



1 



trC 5 



4 
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F16 



with JCjk = 



16 



F17 : 



Edges g*) 



FIB : 



+ ^jib L%n (&jh JCjr* - «4j /Ctm) 

airs of edoes 



paoxs of edqes 



+ L j* L ™ +JC/rn £jbn) 

pairs of edges J 



6 teams. . 
12 terms 

3 terms 



F19 



w 



324 V* 1 V 2 J 



-|2 



+ 



F20 : 

edges 

~ E W*-ft)(^m-&M)JC* m 
pairs of edges 

+ E ( L % ~ (?*) C^™ ~ ^»n) (^/n ^itm + fCj m fC h n ~ Kjh K mn ) 

pairs of edges 
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F21 : 



E Wt-«** 



[ E 9^**) 
E ^ ~ ft) (o^oj - Q/cg - qT^ 4^ J | 



+ 

ftps 



R2 : W 



p = compij, k) 



El : 



m. — -i- - v — - + F. 
' dt " dt 



E2 



MM + CU + KU - R 



E3 
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